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Elastic and inelastic properties of weakly bound s- and p-wave molecules of fermionic atoms that 
collide with a third atom are investigated. Analysis of calculated collisional properties of s-wave 
dimers of fermions in different spin states permit us to compare and highlight the physical mecha- 
nisms that determine the stability of s-wave and p-wave molecules. In contrast to s-wave molecules, 
the collisional properties of p-wave molecules are found to be largely insensitive to variations of the 
p-wave scattering length and that these collisions will usually result in short molecular lifetimes. We 
also discuss the importance of this result for both theories and experiments involving degenerate 
Fermi gases. 
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I. INTRODUCTION 

In the past decade, the advent of magnetic field control 
of interatomic interactions near a Feshbach resonance has 
rapidly sparked both experimental capabilities and the- 
oretical understanding in the field of ultracold quantum 
gases. Using s-wave Feshbach resonances in an atomic 
Fermi gas with two different spin states fundamen- 
tal problems such as the crossover between Bose-Einstein 
condensation and the Bardeen-Cooper-Schrieffer-type su- 
perfluidity [2| have become accessible experimentally, al- 
lowing tests of fundamental theories and studies of novel 
phenomena (see [l| and references within). Recently, 
the observation of p-wave Feshbach resonances in spin- 
polarized Fermi gases [|| instigated theoretical studies 
that predicted a variety of novel many-body phenomena 
Moreover, the recent experimental observation of p- 
wave molecules @ in a 40 K Fermi gas has provided a 
starting point for further investigations of novel many- 
body phenomena in the strongly interacting regime 
The experimental realization of the theoretical predic- 
tions will hinge on the lifetime and stability of p-wave 
molecules, since, e.g., the long lifetime predicted [9j and 
found [l], @ for s-wave molecules was one of the key in- 
gredients that facilitated the observation of fermionic su- 
perfluid behavior in dilute gases. In Ref. [f|, though, it 
was found that p-wave molecules are short-lived which 
limits the possible types of states that can be created for 
p-wave molecular condensates. 

In this paper, we explore the physics of three-body 
collisions involving both s-wave and p-wave molecules 
with the goal of understanding what determines their 
different levels of stability. For both s- and p-wave cases, 
the dominant contribution to the inelastic atom-molecule 
rate coefficient at sufficiently low collision energies has 
s-wave character and therefore does not vanish, even at 
zero temperature. Such processes thus play an important 
role for atom-molecule mixtures near a p-wave Feshbach 
resonance. Our results indicate that many of the charac- 
teristics of p-wave 40 K2 molecules are likely to be shared 
by other atomic species. This expectation is based on 



the universal nature of our results that are discussed be- 
low. We have also found that, near a p-wave Feshbach 
resonance, s-wave atom-molecule collisions are likely to 
be insensitive to the presence of the Feshbach resonance, 
in contrast to the s-wave case. However, our results indi- 
cate that while inelastic p-wave atom-molecule collisions 
are suppressed, elastic p-wave atom-molecule collisions 
are enhanced near a p-wave resonance. We then briefly 
discuss the three-body parameters that might be impor- 
tant for understanding many-body aspects in degenerate 
Fermi gases with resonant p-wave interactions. The pri- 
mary goal of the present study is to discuss calculations 
carried out for both s- and p-wave molecules that high- 
light the main mechanisms that cause their lifetimes to 
be so different. 

Near a two-body Feshbach resonance, the physics that 
determines the collisional properties of s-wave molecules 
is closely related to what has been called Efimov physics 
0, [H ■ Efimov physics occurs when the two-body s-wave 
scattering length, a s = — limfc_>o tan<5 s /fc with 6 S the 
phase shift and k the wave number, is abnormally large 
compared to the characteristic range ro of the interatomic 
interactions. For s-wave molecules composed of fermions 
in different spin states, say FF', Efimov physics predicts 
a universal repulsive effective interaction for FF' + F col- 
lisions which is responsible for the aj 3 ' 33 suppression of 
atomic and molecular losses 8, 9]. For FF' + F collisions, 
Efimov physics also predicts that the atom-molecule scat- 
tering length should be a a d ~ l-2a s 0], which ensures 
a strong atom-molecule interaction in the many-body 
quantum gas, and which provides efficient evaporative 
cooling near a Feshbach resonance. A positive a a d can 
also protect the atom-molecule mixture against collapse. 

However, as has been shown experimentally [5(, near 
a p-wave Feshbach resonance the situation differs strik- 
ingly and p-wave molecules composed of spin-polarized 
fermions in a single internal substate, say FF, tend to 
have short lifetimes. Only recently has atom-molecule 
scattering for large p-wave scattering lengths, a p = 
limfe^ot-tan^p/fc 3 ) 1 / 3 , been studied [H, [l||. (Note 
that the p-wave scattering is sometimes characterized 
equivalently by the "p-wave scattering volume" V p = a 3 
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[LSI Il5l| .) Three-body recombination of spin-stretched 
fermions received earlier attention 0, [l5| where it was 
shown that recombination rate could be significant near 
a p-wave resonance despite the Fermi statistics suppres- 
sion. More recently, the recombination channel in the 
FFF system was studied using a zero-range p-wave pseu- 
dopotential interaction 10]. There, the existence of a 
large number of weakly bound three-body states analo- 
gous to the Efimov states for bosons 0] was proposed. 
In our calculations, we found similar three-body states 
at the simplest level of approximation of the hyperspher- 
ical adiabatic representation. They disappear, however, 
once nonadiabatic effects are included, which leads us 
to speculate that these nonadiabatic corrections might 
also eliminate the three-body states found in Ref. 
Moreover, it is not clear whether the ener gy d ependence 
neglected in their p-wave pseudopotentials fill a ffects the 
universality of the results predicted in Ref. [10] • Our re- 
sults presented here do not suffer from such limitations 
because our two-body finite range model potential intrin- 
sically includes such corrections. 

In this paper we focus on atom-molecule scattering 
properties and find that both elastic and inelastic s-wave 
FF + F processes are largely and surprisingly insensitive 
to variations of a„ and potentially explain the loss rates 
observed in Ref. [g]. This result agrees with the qualita- 
tive arguments in Ref. [l3j. We trace the insensivity of 
the atom-molecule losses to the weakly attractive effec- 
tive atom-molecule interaction and to the fact that the 
size of a p-wave molecule is mainly determined by tq and 
therefore does not change as a p varies, in contrast to s- 
wave molecules whose size is directly related to a s . Our 
results, therefore, indicate that the insensitivity to varia- 
tions of a p of both elastic and inelastic FF + F processes 
is likely to be a universal property of p-wave molecules. 
Based on these results we also speculate on which few- 
body parameters are important for many-body theories 
for spin-polarized Fermi gases. 



II. THEORETICAL BACKGROUND 

We have extracted the three-body elastic and inelas- 
tic collisional properties of s- and p-wave molecules 
from numerically converged solutions of the three-body 
Schrodinger equation carried out in the adiabatic hy- 
perspherical representation [l2| , using model finite-range 
two-body interactions. This representation offers a sim- 
ple, unifying picture from which we can quantitatively 
determine and then qualitatively interpret the origin of 
the significantly different collisional properties of s- and 
p-wave molecules. 

In the adiabatic hyperspherical representation, the 
three-body system is described in terms of the hyper- 
radius R, which gives the overall size of the system, 
and a set of five hyperangles [l2| , which mainly describe 
the interparticle correlations. In this representation the 
Schrodinger equation reduces to a system of coupled or- 



dinary differential equations given (in atomic units) by, 
1 d 2 
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+ V VV ,{R)F V ,{R) = EF U (R), (1) 

where /i = m/y/3 is the three-body reduced mass for 
three identical particles (m being the atomic mass), E 
is the total energy, F v is the hyperradial wave function, 
and v is a collective index that represents all quantum 
numbers necessary to label each channel. In the above ex- 
pression, V vv i are the nonadiabatic couplings calculated 
in terms of the R derivative of the hyperangular solu- 
tions 12]. In Eq. (JTJ) , the nonadiabatic couplings V UV ' 
drive transitions between different channels, represented 
by the effective potentials W v . 

In our calculations we have used the two-body model 
interaction 



v ( r ij) = -£>sech 2 (rv,/r ), 



(2) 



where is the interatomic separation. We have taken 
ro to be the van der Waals length [l6| for 40 K atoms, 
r = ^(mCe) 1 ^ 4 w 65 a.u., and m to be the 40 K atomic 
mass. We varied a s and a p by changing the potential 
depth D around the vicinity of the value that produces 
a weakly bound state of s- and p-wave character, respec- 
tively. For FFF' systems we adjusted D to produce two 
deeply bound FF' molecules (one s-wave and one p-wave) 
in addition to the weakly bound s-wave molecule, and 
we neglected FF interactions. For FFF systems, we ad- 
justed D to produce one deeply bound p-wave molecule 
and one weakly-bound p-wave molecule. In a more re- 
alistic model for p-wave interactions near a Feshbach 
resonance, however, one needs to consider the magnetic 
dipole-dipole interaction which treats the resonance in 
terms of its angular momentum projection m p = and 
Trip — ±1 [171 ]. Our model, which does not include dipole- 
dipole interactions, is applicable to m p — and \m p \ = 1 
resonances individually and does not include coupling be- 
tween m p = and \m p \ = 1 molecules. 

At ultracold energies, the dominant contribution for 
both F + FF and F + FF' collisions comes from the rel- 
ative s-wave, which implies that the dominant three-body 
symmetry for FFF' is = + (J is the total orbital 
angular momentum and ir is the total parity) while for 
FFF it is J 71 " = 1~. Figure [T] shows typical effective 
potentials W v for FFF' [Fig. QTa)] and FFF systems 
[Fig. QIb)] obtained for a s w 300 and a p w 170 a.u., 
respectively. For these scattering lengths, the binding 
energy for s-wave molecules, E§ ~ 1/ma 2 , and p-wave 



molecules, E% 2r p /ma p [18|, are comparable, making 
it easier to contrast the two cases. 

In Fig. Q]the series of effective potentials that converge 
to zero as R — > oo are three-body continuum channels 
representing collisions between three free atoms. Lower- 
lying potentials represent atom-molecule entrance or es- 
cape channels; notice that there are two FF + F channels 



3 




FIG. 1: (Color online) Three-body effective potentials for (a) 
J n — + FFF' fermionic systems with s-wave resonant in- 
teraction and for (b) J w = 1~ FFF spin-polarized fermionic 
systems with p-wave resonant interactions. As is discussed 
in the text, these potentials illustrate the main mechanisms 
that produce different collisional behaviors of the s- and p- 
wave molecules presented here. 



gests that a ad should be proportional to a s , in analogy 
to scattering from a hard sphere of radius a s . The inelas- 
tic processes for this case, however, should be suppressed 
as a s increases, since it becomes increasingly difficult for 
the system to tunnel to the region R ss ro where the cou- 
pling to the deeper molecular states lies [we have verified 
this statement from our calculations of the nonadiabatic 
couplings in Eq. |(TJ)]. fn fact, in Ref. [8] it was demon- 
strated that the suppression of the inelastic losses scale 
with a as a]~ 2pa — a,T 3 - 33 , where po is related with the 
strength of the repulsive potential barrier in the range 
ro <C R <C a s - In contrast, the presence of an attrac- 
tive potential for FF + F collisions instead of a repul- 
sive potential allows the atom and molecule to approach 
each other closely without any suppression of the inelas- 
tic transition probability. In addition, the absence of a 
repulsive barrier does not allow us to make an analogy 
to scattering from a hard sphere, which implies that a ad 
is not necessarily a simple function proportional to a p . 

In Refs. [HI, it has been speculated that there 
is one trimer state in this attractive potential that can 
cause resonant enhancement of the elastic and inelastic 
rates when varying a p . As we show below, in our cal- 
culations (see Fig. we do observe resonant effects for 
elastic and inelastic FF + F collisions, and we do asso- 
ciate them with the formation of a trimer state as pre- 
dicted in Refs. [H,[l3|- However, we believe such trimer 
states are likely to be non-universal in the sense that 
their energy as a function of a p depends on the short- 
range physics not fully represented by the two-body pa- 
rameters in Ref. fl3l Ibj]. Therefore, the resonant peaks 
in the elastic and inelastic rates will likely depend on 
a purely three-body non-universal parameter, similar to 
the problem with identical bosons where all three-body 
observables depend on a three-bo dy p arameter represent- 
ing the short-range physics fft, M, fl9f. 



III. ELASTIC AND INELASTIC 
ATOM-MOLECULE COLLISIONS 



in Fig.[TJb) converging to the p-wave FF molecule energy 
at large R due the possible I — and I = 2 atom-molecule 
relative angular momenta allowed for J n = 1~. In Fig. [1] 
we also sketch the elastic and inelastic pathways that are 
of interest here. 

From the figure, we notice a fundamental difference be- 
tween the initial FF' + F and FF + F collision channels. 
For FF' + F collisions, the effective potential is repulsive 
(proportional to 1/R 2 @) in the range ro <C R <C a s , 
while the I = effective potential for FF + F colli- 
sions is attractive in the analogous range r <C R <C 

a r>^ 2 / r p^ 2 [the quantity a p ^ 2 /r^ 2 can be understood from 
the length scale denned by the molecular binding energy 
l/k p — (m£j) -1 / 2 ]. The presence of a repulsive barrier 
for FF' + F collisions whose range scales with a s sug- 



We have verified the qualitative behavior described 
above for the elastic and inelastic rates by solving the 
Schrodinger equation |1| in coupled hyperradial form us- 
ing an R- matrix approach [20|. We define the atom- 
molecule scattering length in terms of the real part of the 
complex phase-shift obtained from the S-matrix element 
[exp(2i6 ac i) — S a d,ad] associated with the atom-molecule 
channel, 

Re[tan<5 ad ] 
aw = - km , (3) 

kad^O k ad 

with k ad = [2fi ad (E + Eb)} 1 / 2 where fi ad = 2to/3 and E b 
is the binding energy of the weakly bound dimer. For 
inelastic collisions, we define the vibrational relaxation 
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rate V Ie i as 



V„ 



E 



(2J+ 1)tt 



(4) 



where i represents the initial collision channel 
(atom+weakly bound molecule); and /, all possible 
final collision channels (atom+deeply bound molecule). 
We have calculated a a d and V rc \ for values of a s and a p 
ranging from r (~ 65 a.u.) up to approximately 70ro at 
a collision energy of 1 nK, in order to satisfy the k ac i — > 
limit in Eq. ([3]). The results are shown in Figs.[2ja) and 
(b). 

For FF' + F collisions, we have obtained both a a d and 
V^ei in the a s r regime that verifies the a a d ~ 1.2a s 
universal prediction [9( and reproduces the strong a,T 3 ' 33 
suppression of the inelastic rate @, Q . This scenario is 
rather favorable for the realization of cold atom-molecule 
mixtures since it combines strong elastic rates with sup- 
pression of atomic and molecular losses. 
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FIG. 2: (Color online) Numerical calculations for the (a) 
atom-molecule scattering length a a d and (b) vibrational re- 
laxation rate V Ic i for J n = + FF' + F (dash-circle line) and 
J n — l~ FF + F (solid-circle line) collisions as a function of 
as and a p , respectively. While for FF' + F collisions a a d and 
V le \ show a strong dependence on a 3 , for FF + F they are 
essentially insensitive to variations of a p . 



On the other hand, FF + F elastic and inelastic pro- 
cesses involving p-wave molecules are quite different. For 
smaller values of a p , we observe a strong variation of both 
a a d and V re i [see Fig. [5] due to the presence of a trimer 
state [ill, that becomes unbound as a p increases. As 
we mentioned before, we do not expect the energy of such 
a trimer state, nor the value of a p at which it appears, 
to be universal, since they both depend on the short- 
range physics. Consequently, the resonance position in 
Fig. can change, and eventually disappear, for differ- 
ent choices of the two-body interaction. However, near a 
p-wave two-body Feshbach resonance, i.e., in the regime 
Op 3> ro, our results for a a d and V le i are insensitive to 
variations of a p . We rationalize this result by recogniz- 
ing that even though a p changes substantially, the size 
of the weakly-bound p-wave molecule is practically un- 
changed. We expect that this insensitivity of a a d and 
V re i to variations of a p is universal but that their actual 
values will not be universal since the system must access 
small R where details of the two-body interactions be- 
come important. Further, it can be demonstrated by a 
simple WKB analysis [|[ that the constant character of 
V Te \ as a function of a p is ensured provided the effective 
potential falls off faster than 1/R 2 . From our numeri- 
cal calculations, for instance, we found that a a d < in 
the regime a p 3> ro, in contrast to the a a d > result of 
Refs. Ill, [3. 

The consequences of the three-body physics discussed 
above for the lifetime and stability of spin-polarized 
fcrmionic mixtures of atoms and molecules are very dif- 
ferent from the scenario found for a gas of fermions in 
two different spin states. The constant value of a a d for 
o-p ^ fo means that atom-molecule collisions will not 
help much for evaporative cooling, unless a a d is unnat- 
urally large. More importantly, the insensitivity of V rc \ 
to a p implies that the molecular lifetimes do not change 
as we approach the point of divergence of a p , in con- 
trast to s-wave molecules whose lifetime increases near a 
zero-energy Fano- Feshbach resonance. These results are 
in agreement with recent experimental data for an atom- 
molecule gas mixture of 40 K atoms Q , where it was found 
that the lifetimes of both m p = and \m p \ = 1 molecules 
do not depend on the magnetic field and therefore on a p . 
In Ref. [H , however, it was found that the molecules were 
shorter-lived than one would expect based on the mag- 
netic field independent two-body dipolar relaxation rates. 
More recently [2l[, loss rates substantially below those 
found in [j| were observed after the atoms had been re- 
moved from the trap, leaving only molecules. Therefore, 
since both dipolar relaxation and atom-molecule colli- 
sions yield molecular losses independent of the magnetic- 
field, we believe that the shorter molecular lifetimes were 
due to vibrational relaxation as we have calculated here. 
Our prediction of approximately constant atom-molecule 
losses also demonstrates that even when atoms are pre- 
pared in their lowest hyperfine state, and therefore dipo- 
lar relaxation is energetically forbidden, a magnetic field 
independent molecular lifetime is likely to be found for 
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other atomic species, but now limited by three-body loss 
processes. In the absence of atoms, molecule-molecule 
elastic and inelastic processes should also be important. 
We speculate that these are also likely to be insensitive 
to variations of a p (see also Ref. (lij). based on the fact 
that the size of a p-wave molecule does not depend on 
a p . 

IV. FEW-BODY PARAMETERS FOR 
MANY-BODY THEORIES 

The fact that a a d does not seem to depend on a p raises 
the question whether a a d will be the important param- 
eter for theories attempting to describe many-body as- 
pects of atom-molecule mixtures in spin-polarized Fermi 
gases. A constant value of a a d across a Feshbach res- 
onance essentially means that the atom-molecule inter- 
action does not "feel" the presence of the p-wave reso- 
nance. From this perspective, the atom-molecule inter- 
action is not controllable, and the atom-molecule mix- 
ture cannot be claimed to be strongly interacting near 
the resonance. For this reason, we have also performed 
some preliminary calculations to explore p-wave atom- 
molecule collisions. For J 71 " = + , 1 + , and 2+ (which are 
all the symmetries that allow p-wave collisions between 
an atom and a p-wave molecule) , we have found that the 
three-body effective potentials are now repulsive in the 
range tq <C R <C a p ^ 2 /r]/ 2 , instead of attractive as found 
for J 71 " = 1~ . The presence of such a repulsive barrier 
suppresses inelastic p-wave atom- molecule collisions Q, 
making these collisions less important than s-wave col- 
lisions in determining the lifetime of p-wave molecules. 
On the other hand, the repulsive barrier indicates that 
the p-wave atom-molecule scattering length, a p ad , should 

be proportional to ap^ 2 /r^/ 2 , again in analogy to scat- 
tering from a hard-sphere. Therefore, it is likely that 
a p ad will be an important parameter for many-body theo- 
ries of atom-molecule mixtures near a p-wave resonance, 
since aP ad is a controllable parameter. Experimentally, 
however, p-wave atom-molecule interactions will be com- 
parably difficult to control as p-wave atom-atom interac- 
tions because of the suppression of the energy dependent 
elastic cross-sections. Certainly, it will be necessary to 
explore in depth the importance of a v ad , but this discus- 
sion is outside the scope of this paper. 



V. SUMMARY 



In this paper, we have determined the atom- molecule 
scattering length, a a d, and vibrational relaxation rate, 
Vrci, for collisions involving s- and p-wave molecules of 
fermionic atoms as a function of a s and a p which can be 
controlled with two-body Feshbach resonances. We have 
determined a a d and V le \ for s-wave molecules of fermionic 
atoms in different spin states, verifying previously de- 
rived results. We could thus identify the mechanisms that 
control the collisional properties of p-wave molecules. For 
spin-polarized systems, we have found that both a a d and 
V Te \ are insensitive to variations of the two-body p-wave 
scattering length a p for a p 3> rn, consistent with the re- 
cent experimental data in Ref. [f| . We expect the results 
for a a d and V Ye i to be universal with respect to their in- 
sensitivity to variations of a p but their values can change 
for different atomic species. We also speculated that, due 
to the constant character of a a d, the relevant parameter 
for many-body theories for a spin-polarized mixture of 
atoms and molecules is likely to be associated with p- 
wave collisions between a fermionic atom and a p-wave 
molecule. Our preliminary results indicate that the p- 
wave scattering length for such collisions, a p ad , is likely to 

be proportional to a p ^ 2 /r p ^ 2 . 

Note that the very recent preprint by Jona-Lasinio et 
al. [HI argues that the threshold exponent for three- 
body recombination of spin-polarized fermions should be 
proportional to K3 oc a p 5 ^ 2 r p ^ 2 at large p-wave scattering 
lengths. This differs slightly from the prediction of Suno 
et al. ,15] that K 3 oc a*. The derivation of Suno et 
al. was, in fact, based on dimensional analysis, starting 
from the assumption that a p was the dominant length 
scale. In fact, it is now understood that a p alone does not 
adequately describe the dominant length scale for two- 
body p-wave scattering, implying that a modification of 
the dimensional analysis prediction of Suno et al. should 
be expected. Investigation of this point will be relegated 
to future studies. 
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